New constraints on neutrino physics from Boomerang data 
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Steen Hannestad 
NORDITA, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
(May 1, 2000) 

We have performed a likelihood analysis of the recent data on the Cosmic Microwave Background 
Radiation (CMBR) anisotropy taken by the Boomerang experiment. We find that this data places 
a strong upper bound on the radiation density present at recombination. Expressed in terms of the 
equivalent number of neutrino species the 2a bound is N v < 13, and the standard model prediction, 
N v — 3.04, is completely consistent the the data. This bound is complementary to the one found 
from Big Bang nucleosynthesis considerations in that it applies to any type of radiation, i.e. it is 
not flavour sensitive. It also applies to the universe at a much later epoch, and as such places severe 
limits on scenarios with decaying neutrinos. The bound also yields a firm upper limit on the lepton 
asymmetry in the universe. 

PACS numbers: 98.70.Vc, 14.60.St, 13.35.Hb 



Introduction — The standard Big Bang model has been 
remarkably successful in describing the observed features 
in our universe [jjj. The latest and most impressive 
confirmation of the model comes from observations of 
anisotropics in the cosmic microwave background radi- 
ation (CMBR). These anisotropics are predicted by the 
Big Bang model, and were first discovered by the COBE 
satellite in 1992 Subsequently it was realized that 
precision measurements of these anisotropics can yield 
very detailed information about the fundamental cosmo- 
logical parameters |^] , and accordingly a vast number of 
papers have investigated the potential of upcoming ex- 
periments to measure these parameters (see for instance 
and references therein). 

Now we have the first results which may rightly be 
called precision CMBR measurements. They stem from 
the balloon-borne experiment Boomerang which was 
flown over the Antarctica in 1999 ||. 

The results indicate that the universe is flat, and are es- 
sentially a confirmation of the standard Big Bang model 

The data is so good that it can also be used to constrain 
physics beyond the particle physics standard model. A 
strong indication of such exotic physics would be addi- 
tional radiation energy in the universe at the time of re- 
combination. This could for instance be caused by addi- 
tional light neutrinos, or some other exotic particle which 
decoupled at very high temperature (jjj. In the present 
letter we use the data from Boomerang to place a strict 
upper limit on the radiation density present at the time of 
recombination (T ~ 1 eV) . The standard way of express- 
ing the energy density in light, non-interacting species, 
is in terms of the equivalent number of neutrinos |IJ 



N eB}i = 



Pi 



(1) 



where p„.o is the energy density in a standard neutrino 
species. 



From Big Bang nucleosynthesis one can also infer a 
limit to the effective number of neutrino species. By ob- 
serving the primordial abundances of D, 4 He and 7 Li, 
and comparing them to the theoretically predicted val- 
ues, one can infer an upper limit to AT„ of Jl^] 

A^.bbn < 4. (2) 

The CMBR limit can be viewed as complementary to 
the BBN limit because the limit from BBN applies to the 
radiation energy density present when the temperature of 
the universe was about 1 MeV, whereas the CMBR limit 
applies at a temperature of 1 eV. Furthermore, the BBN 
limit is also flavour sensitive. If the extra energy density 
is in the form of electron neutrinos, it changes the weak 
reaction rates for the beta processes that interconvert 
neutrons and protons. With some fine tuning, even a 
very substantial amount of energy can be hidden in the 
neutrino sector while still yielding the same outcome for 
BBN. 

That is not the case for CMBR. In this case, extra en- 
ergy density is detectable because the CMBR spectrum 
changes with the addition of radiation. After recombi- 
nation, the CMBR fluctuations can still change. If the 
universe is completely matter dominated and flat, the 
photons see a constant gravitational potential (in the lin- 
ear approximation), and thus travel with constant energy, 
except for the overall redshifting. However, immediately 
after recombination the universe was not completely mat- 
ter dominated so that the gravitational potential was not 
constant. This leads to an enhancement of the first acous- 
tic peak in the power spectrum and is known as the early 
Integrated Sachs- Wolfe (ISW) effect §. This effect is 
only sensitive to energy density and not to the specific 
nature of the radiation. 

Likelihood analysis — The CMBR fluctuations are usu- 
ally described in terms of the power spectrum, which is 
again expressed in terms of C; coefficients as 1(1 + 1)C/, 
where 



1 



Ci = (\a lm \ 2 ). 



(3) 



The ai m coefficients are given in terms of the actual tem- 
perature fluctuations as 



T(6,<f>) = Y,a lm Y lm {6,, 



(4) 



Given a set of experimental measurements, the likelihood 
function is 



£(6) ocexp ( -^[Ctey^x 



(5) 



where = (f2, SI;,, Hq, n, r, . . .) is a vector describing the 
given point in parameter space, a; is a vector containing 
all the data points and C(9) is the data covariance ma- 
trix. This applies when the errors are Gaussian. If we 
also assume that the errors are uncorrelated, this can be 
reduced to the simple expression, C oc e~ x ' , where 



(C(,obs — Cl t theory)i 



(6) 



is a x 2 -statistics and iV max is the number of power spec- 
trum data points |]l3| . In the present letter we use Eq. (||) 
for calculating x 2 - 

The procedure is then to calculate the likelihood func- 
tion over the space of cosmological parameters. The like- 
lihood function for N v is then obtained by keeping N v 
fixed and maximizing £ over the remaining parameter 
space. The fundamental free parameters which we al- 
low to vary are: Jl m , the matter density, fib, the baryon 
density, h, the Hubble parameter, n, the spectral index, 
r, the optical depth to reionization, and Q, the over- 
all normalization of the spectrum given in terms of the 
quadrupole moment 0. The range in which they are al- 
lowed to vary are given in Table I. We assume a fiat 
universe with Oo = fit + f^m + = 1- This is the 
value strongly suggested by the Boomerang experiment 
Also, the final result of the analysis does not vary 
much even if Jlo is allowed to vary. For a given value of N v 
we maximize the likelihood over the remaining parame- 
ter space by using the non-linear optimization method 
called simulated annealing jL4|. The data set which we 
use is the one given in de Bernardis et al. || , and we use 
the publicly available CMBFAST package for calculating 
theoretical power spectra ]l5| . 

We show the result of the likelihood analysis in Fig. 1 
in terms of % 2 . The absolute minimum \ 2 is 5.27 at 



*Note that there is an estimated 10% calibration uncertainty 
in the overall normalization of the Boomerang data. However, 
this effect is completely degenerate with varying Q because 
we use only one data set. Therefore we do not need to worry 
about it. 



TABLE I. The free parameters used in the present analysis, 
as well as the range in which they are allowed to vary. 



Parameter 



range 



Q 

Q b h 2 

h 

n 



5-40 nK 
0.2-1 
0.002-0.050 
0.30-0.9 
0.7-1.3 
0-1 



N v = 1.6. The number of Boomerang data points is 
12 and we allow 6 parameters to vary. Thus, the 
number of degrees of freedom in our fit is 6, so that the 
obtained best fit x 2 falls within la of the expected for an 
acceptable fit. The CMBR constraint on N v is given by 



N„ < { 



6.2 
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la, 
2a. 



(7) 



While this constraint is clearly much weaker than the 
N p < 4 obtained from BBN considerations |]l2| , it applies 
to any type of relativistic energy density. 




FIG. 1. x f° r the Boomerang data as a function of N v . 
The curve has been obtained by minimizing \ 2 over all other 
free parameters. 

Discussion — We have calculated a strong upper bound 
on radiation density present at recombination. Next, we 
may ask what bounds can be placed on neutrino proper- 
ties from this. 

The standard model prediction is that N u = 3.04 Jl6| . 
Although the absolute minimum for \ 2 is at N v = 1.6, 
the standard model is completely consistent with the 



2 



Boomerang data. There is no indication in the data of 
neutrino physics beyond the standard model. 

If sterile neutrino degrees of freedom exist, then they 
can be exited by oscillations in the early universe 0,|l8| • 
However, a sterile neutrino can at most contribute 
■^sterile — 1 so that our CMBR bound does not even 
exclude this possibility at the la level. 

Next, we have no direct measurement of the lepton 
asymmetry in the universe, and quite a large lepton 
asymmetry could in fact be hidden in the neutrino sector 
[ pToj p2| . The lepton asymmetry in neutrinos is usually ex- 
pressed in units of £„ = \i v jT v . The neutrino distribution 
function is then given by / = (exp(E/T v ± £„) + 1) — 1 , 
where + applies to neutrinos and — to antineutrinos. 
From BBN, one obtains the bound [EOl 



6 [-0.06,1.1] 
|£wl<6.9. 



(8) 
(9) 



This is in the absence of any oscillations. For massless 
neutrinos, a chemical potential is equivalent to a change 
in the effective number of species [GlJ 
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Using this, we can translate our 2a upper bound on N v 
to a bound on £„ 



< 3.7 (2a). 



(11) 



This bound applies if only one species has a chemical 
potential. If the asymmetry is equally shared then |£„| < 
2.4. Lesgourgues and Peloso Q have also discussed a 
cosmological lepton asymmetry as a possible explanation 
of the relatively low amplitude of the second acoustic 
peak in the Boomerang data, but they did not perform 
a likelihood analysis of the data. 

Interestingly this bound is at the brink of excluding a 
scenario for the production of ultra high energy cosmic 
rays (UHECRs), proposed by Gelmini and Kusenko 
In this scenario, very high energy neutrinos interact with 
a degenerate sea of neutrinos within about 50 Mpc of 
earth. These interactions produce charged particles such 
as protons which would then be the observed UHECR 
primaries. However, because of the small vv interaction 
cross section this scenario only works if the cosmic neu- 
trino background is degenerate with £ ~ 4. At present, 
the CMBR bound disfavours this model, but cannot ex- 
clude it completely. It has been estimated that with the 
upcoming MAP and Planck experiments, it will be possi- 
ble to constrain £ with a precision of about 0.1 j22|. This 
will definitively confirm or rule out scenarios like this. 

Neutrino decays to massless secondaries prior to re- 
combination are also seen in the CMBR fluctuations as 
an increased N v p5]-[29|] . The effective number of neutri- 
nos is given roughly by 



N u a 3 + 0.516a 2/3 , (12) 
where a is the decay "relativity" parameter 

Again we can translate our bound on N v into a bound 
on a and thereby on neutrino decays 



a < 85.3 (2a). 



(14) 



In Fig. 2 we show what this translates into in terms of 
neutrino lifetime and mass. 
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FIG. 2. The excluded region for neutrino decays. The ho- 
risontal line indicates the time of recombination. Decays tak- 
ing place after this are not constrained by the present calcu- 
lation. The other full line is from Eq. (|l4|). 

Note that CMBR measurements can also be used to 
constrain late neutrino decays, which take place after 
recombination p6|-p9|. However, such decays are de- 
tectable because of the late ISW effect they produce at 
small I. Since Boomerang does not detect fluctuations 
below about I — 50, we have not calculated any bound 
on late decays. 

In conclusion, Boomerang has provided us with the 
first precision CMBR data. This data can be used to 
constrain many different types of exotic physics, not just 
in the neutrino sector. The Boomerang data is only the 
first indication of what we can expect with the new gen- 
eration of CMBR experiments. In the next few years we 
will have data from the MAP and Planck satellites |30|, 
which is expected to be an order of magnitude better 
than that from Boomerang. 
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